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We derive a rigorous theory of the interaction between photons and spatially extended excitons 
confined in quantum dots in inhomogeneous photonic materials. We show that, beyond the dipole 
approximation, the radiative decay rate is proportional to a non-local interaction function, which 
describes the interaction between light and spatially extended excitons. In this regime, light and 
matter degrees of freedom cannot be separated and a complex interplay between the nanostruc- 
tured optical environment and the exciton envelope function emerges. We illustrate this by specific 
examples and derive a series of important analytical relations, which are useful for applying the 
formalism to practical problems. In the dipole limit, the decay rate is proportional to the projected 
local density of optical states and we obtain the strong and weak confinement regimes as special 
cases. 
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I. INTRODUCTION 

The dipole approximation (DA) is one of the most cen- 
tral and successful approximations in quantum optics and 
quantum electrodynamics (QED). When describing the 
Hght-mattcr interaction, the DA is valid if the variation 
of the electromagnetic field is negligible over the spa- 
tial extent of the emitter. Since optical wavelengths ex- 
ceed atomic dimensions by orders of magnitude, this is 
an excellent approximation in atomic physics. The ad- 
vances in solid-state quantum optics have enabled the re- 
alization of semiconductor nanostructures with strongly 
modified optical properties and embedded self-assembled 
quantum dots (QDs). Both atoms and QDs have a dis- 
crete spectrum with optically active transitions but as 
opposed to atoms, QDs are inherently mesoscopic solid- 
state structures whose transition energy, position, and 
chemical composition can be controlled by semiconduc- 
tor nanotechnology, and in nanophotonic structures the 
electromagnetic environment can have pronounced spa- 
tial variations. For QDs the validity of the DA is not clear 
a priori, and the purpose of this article is to derive the 
theory of spontaneous emission beyond the DA for exci- 
tons confined in QDs embedded in nanostructures. Re- 
cently, a large deviation from dipole theory was observed 
for small QDs in close proximity to a metallic mirror—, 
directly illustrating the need for a theory of light-matter 
interaction beyond the DA. 

The tunability of QD sizes can lead to very interesting 
exciton effects. For strongly confined states in small QDs 
the Coulomb interaction can be neglected and excitons 
can thus be described as mutually independent electrons 
and holes^i^. For larger QDs the Coulomb interaction 
plays an increasingly important role and the electron and 
the hole form a bound exciton state, which has an oscilla- 
tor strength (OS) proportional to the volume of the exci- 



ton. This is the so-called giant-OS effect^"™, which has re- 
ceived particular attention in the field of solid state QED 
because of the theoretical prediction^, that QDs must be 
in the giant- OS regime in order to achieve strong coupling 
between a single QD and a microcavity. Indeed, some of 
the first demonstrations of strong coupling in microcavi- 
ties were achieved with large QDs^i^. However, the rapid 
increase in the quality factors and the reduction of mode 
volumes in photonic crystal membrane nanocavities over 

the past years have enabled strong coupling using small 
QD9ia-i2. 

A key signature of the giant-OS effect is fast radia- 
tive decay ratesi^. Fast total decay rates have been 
observedi^^ but the non-radiative decay rate was not 
measured in any studies of large QDs except for a re- 
cent worki^. There it was shown that contrary to the 
common assumption, non-radiative recombination can be 
the dominant decay process for large QDs resulting in a 
small OS but a fast total decay rate. Therefore, measur- 
ing the non-radiative decay rate is essential to drawing 
conclusions about the OS and this has so far hindered 
a complete experimental demonstration of the giant-OS 
effect. 

Another effect of increasing the QD size, which has 
received much less attentioni^"— , is the fact that for suf- 
ficently large QDs the DA may break down. The usual 
criterion for the validity of the DA is that the product 
of the length of the optical wave vector k and the spa- 
tial extent of the emitter £ must be much smaller than 
unityi^, i.e., |k|£ ^ 1, but this criterion is insufficient 
to ensure the validity of the DA for QDs in nanostruc- 
tures. For QDs, there are four reasons why the DA could 
break down. Firstly, the QDs are embedded in semi- 
conductors with a high refractive index, e.g., n « 3.4 
for GaAs, which increases |k|. Secondly, QDs can be as 
large as £ = 100 nm in lateral size^. As an example, for 
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FIG. 1. Comparison between the LDOS and the non-local in- 
teraction function introduced in this work. An exciton emit- 
ter with a given envelope function (dark gray) is embedded in 
an inhomogeneous dielectric environment indicated schemat- 
ically by the different refractive indices ni (white) and 712 
(light gray). The Green's tensor G{r,r',uj) is a propagator 
of the electric field between two spatial points at a given fre- 
quency uj and it is depicted as the arrows. (A) In the DA the 
spontaneous emission rate is governed by the LDOS, which is 
given by the imaginary part of the Green's tensor evaluated at 
(ro, ro), where ro is the center of the emitter. (B) Beyond the 
DA the spontaneous emission rate is governed by a non-local 
interaction function, which is given by an integral over the 
imaginary part of the Green's tensor connecting all possible 
combinations of r and r', weighted by the envelope function. 



£ = 45 nm and a free-space wavelength of 970 nm we ob- 
tain |k|£ « 1. Thirdly, the criterion stated above is valid 
for homogeneous media. In nanostructures the optical 
field modes can have strong gradients rendering the DA 
invalid even if |k|£ ^ 1 is fulfilled in bulk for light at 
the same frequency. Fourthly, the influence of the finite 
size of the QDs is enhanced by the asymmetric nature of 
exciton wave functions in QDsi. Thus, a proper theory 
of spontaneous emission beyond the DA must be valid 
for arbitrary electromagnetic environments. Here we de- 
rive such a theory from first principles and show that the 
radiative decay rate depends on a non-local interaction 
function, whose physical interpretation is illustrated in 
Fig. [TJ In the DA (Fig. [HA)) the radiative decay rate 
is proportional to the projected local density of optical 
states (LDOS), which describes the field amplitude at 
the position of the emitter due to emission from a dipolc 
source at the same position. Beyond the DA (Fig.[IjB)), 
a double integral over all points in space must be per- 
formed, where the integrand is weighted by the envelope 
function of the exciton. The double integral describes 
the physics emerging beyond the DA. These effects find 
a natural description within the framework of the elec- 
tromagnetic Green's tensor. Thus, the self-interference 
giving rise to spontaneous emission is described mathe- 
matically by the imaginary part of the electromagnetic 
Green's tensor. 

The non-local aspect of light-matter interaction be- 
yond the DA implies that light and matter degrees of 
freedom cannot be separated, i.e., the radiative decay 
rate is neither proportional to the projected LDOS nor 
to the OS. This points to another previously overlooked 
problem in the interpretation of the experimental results 
on large QDs: even if the non-radiative decay rate dis- 



cussed above had been measured and found negligible, 
the highly non-trivial influence of the ubiquitous sur- 
rounding optical nanostructure cannot be approximated 
by a homogeneous medium and the OS has no general 
physical meaning. Here we present the complete quan- 
tum theory of spontaneous emission for two-level QDs 
in inhomogeneous media, which provides the theoretical 
framework for more quantitative future experiments and 
enables calculating non-Markovian decay dynamics and 
radiative (Lamb) shifts. We consider InGaAs QDs, but 
our formalism can be readily modified to describe other 
materials. 

This paper is organized as follows: In section I we de- 
scribe the exciton state and in section II we calculate 
the spontaneous emission from cxcitons beyond the DA 
in the Wigner-Weisskopf model. In section HI we de- 
rive the connection to dyadic Green's tensors, introduce 
a non-local interaction function and discuss the physical 
implications of the results. We consider the dipole limit 
of our results in section IV and we apply our formalism 
to three special cases in section V. Finally, we present the 
conclusions in section VI. In Appendix A we solve the ef- 
fective mass equation for the geometries relevant for this 
work. In Appendix B we consider the classical analogue 
of our results. The derivation of the relation to Green's 
tensors is included in Appendix C. Finally, in Appendix 
D we show the analytical calculation of the decay rate of 
spherical excitons beyond the DA. 



II. EXCITONS IN QUANTUM DOTS 



A bulk semiconductor consists of nuclei and electrons 
and in the Born-Oppenheimer approximation the motion 
of nuclei and electrons are decoupled. Thus we consider 
electron states imposed on an equilibrium state of the 
nuclei. At low temperatures the valence bands are com- 
pletely filled while the conduction bands are empty. Al- 
though the upper valence bands in InAs and GaAs are 
degenerate in bulk materials, this degeneracy is lifted in 
the presence of confinement and strain in QDs and we 
consider the heavy-hole band only. This is a good ap- 
proximation for low excitation powers and low tempera- 
turesSS. 

The description of the ground state of a bulk semi- 
conductor consisting of N electrons in the upper valence 
band must be treated in a many-body formalisro^i^^— . 
In the simplest possible case, i.e., when neglecting inter- 
actions, the ground state wave function is given by the 
Slater determinant^. We neglect the spin degree of free- 
dom, which amounts to considering only bright excitons 
in which the electron and hole spins are antiparalle l^^'^^ . 
It is convenient to use the compact occupation number 
formalism and we define the ground state of the crystal 
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as the Fermi sea \T) given by 



|lv,ki, • ■ 



1 lv,kiv) 



(1) 



where Cv,ki {c\ ^ ) annihilation (creation) operator 

of an electron in the valence band with k = k,; and \Ojr) 
denotes the state void of any electrons. These second- 
quantization operators create single-particle states with 
corresponding wave functions 'i/'v.kilrOj which may be 
written in Bloch-form as 



V 



e''''''''Uv,k,(ri 



(2) 



Here the valence band Bloch function Wv.ki(rj) has the 
periodicity of the crystal lattice and is normalized over a 
unit cell and V denotes the crystal volume. 

We can write an excited state of the bulk semiconduc- 
tor as 



kckv 



'^c,k 



Cv.kvl-^) 



(3) 



where Cc,kc (c^ k ) annihilation (creation) opera- 

tor of an electron in the conduction band with k = kc. 
The operators considered above are in the electron rep- 
resentation, but at this point it is convenient to change 
to the electron-hole representation by defining the fol- 
lowing operators^! = Cc,k, 



,kv' 



and 6]^^ = Cv,kvj where the electron operators have sim- 
ply been renamed and 6kh (^kh) denotes the annihilation 
(creation) operator of a hole in the valence band. With 
this convention we write the excited state of the bulk 
semiconductor as 



(4) 



This definition of the electron-hole representation has a 
number of consequences for the properties of holes. In 
particular, the following transformations hold, where the 
subscript h refers to holes in the electron-hole representa- 
tion and the subscript v refers to electrons in the valence 
band in the electron representationSiiSS: k^ = — kv (wave 
vector), -Bh.v = ~Ec_v (energy), uih = —m^ (effective 
mass), (7h = —q (charge), and Vh(r) ~ —Vv{r) (confine- 
ment potential). Here i?h.v (Ec.v) denotes the energy of 
a hole (an electron) relative to the valence band edge en- 
ergy Ev and q is the negative of the elementary charge, 
i.e., q = -\e\. 

In the presence of Coulomb interaction and/or quan- 
tum confinement potentials the states |Xk„kh) are no 
longer eigenstates. Instead the new exciton eigenstate 
\X) can be expanded as 



l-^> = Xke,kh|Xkekh) 

ko.kh 



(5) 



where Xk^Mh are expansion coefficients. The correspond- 
ing wave function X(ro,ro,rh) can be found by projec- 
tion onto the position eigenvectors, where we have ex- 
plicitly included the center position of the QD, Vq, 

X(ro,rc,rh)= ^ Xk<,,kh (rcrhl^kekj 



kc.kh 



Xko.kh 



ko.kh 



X e 



ikcTc 



,k,(re)e"^-"-''iiv,k,(rh) (6) 



-X(i"o, Tc, rh)wc,o(i'c)uv,o(i"h) 



(7) 



where Bloch's theorem, the transformation of a sum to 
an integral, — >■ j^jw J dk, and the definition of an in- 
verse Fourier transform have been used. The last equal- 
ity holds when only excitations near the band edge are 
considered so that the Bloch functions may be evalu- 
ated at k = 0, which is a good approximation for low 
temperatures and low excitation intensities. The func- 
tion x(ro,re,rh) is denoted the exciton envelope func- 
tion, which is given by the solution to the effective mass 
equation2i 

i/EM(ro,rc,rh)x(ro,rc,rh) = (E - Eg)x{ro,r^,rY,), (8) 

where E is the exciton energy. Eg = Ec — E^ is the band 
gap energy, and the effective-mass Hamiltonian is given 
by 



i7EM(ro,ro,rh) 



-F Vh(ro,rh) - 



K=(ro,rc) 



q 



(9) 



47reoei.|ro - rh| 



Here Po (ph) is the electron (hole) momentum operator. 
Too is the electron rest mass, eg is the vacuum permit- 
tivity, and Ej- is the relative dielectric constant. In the 
strong-confinement limit where the Coulomb interaction 
can be neglected the solution to Eq. is x(ro,rc,rh) = 
^;(ro,ro)-Fh(ro,rh), where Fc(ro,rc) and Fh(ro,rh) de- 
note the electron and hole envelope functions, respec- 
tively. In the weak- confinement limit, electrons and holes 
are entangled and therefore their wave functions do not 
separate. We consider solutions to the effective-mass 
equation in specific geometries in Appendix [X] 



III. QUANTUM THEORY OF SPONTANEOUS 
EMISSION BEYOND THE DIPOLE 
APPROXIMATION 

We describe light-matter interaction by the minimal 
coupling Hamiltonian in the generalized Coulomb gauge 
in which we assume V • (er(r)A(r, t)) = 0, where A(r,i) 
is the vector potential. The interaction Hamiltonian 
reads22*30 



i/'(r,t)=^A(r,0-V. 

Too 



(10) 



4 



The vector potential is given h}dl. 



A(r,i) 



E 



(Mr) 



(11) 



where fj, = (k, s) is the combined wavevector k and po- 
larization index s € {1, 2}, oj^ is the optical angular fre- 

V ^ 



quency, e. 



is a normalization constant, is 

the polarization unit vector, and ^^i(r) is the field distri- 
bution function that solves the vector Helmholtz equa- 
tion with fixed boundary conditions, and aj^ are the 
field annihilation and creation operators, respectively. In 
second quantization the interaction Hamiltonian can be 
written as 



H' 



a,l3 



(12) 



-iujct 



where «,/? G {c,v} and H'^^^-^, = (l^^k, |^^'(r, i)|l/3,k,)- 

It is convenient to use the interaction picture where 
the time-evolution of the operators is governed by the 
non-interacting Hamiltonian. When considering only ex- 
citations near the band edge k = 0, the energies of the 
conduction and valence bands and hence the transition 
energy of the QD, hiOQ, do not depend on i or j. The oper 
ators in the interaction picture become c^^k = Cck^ 
elk = cLe'"^-^*, Cv,k = Cv,ke-''^-*, and k 
When inserting these operators in Eq. we obtain 
terms proportional to e±i('^f<+'^o)t ^j^jj ^±iu^t^ which are 
rapidly oscillating as a function of time as well as the 
slowly oscillating terms proportional to e*"^''* where 
~ t^/i — t^o- The rapidly oscillating terms average 
to zero and are therefore neglected in the rotating wave 
approximation. 

We consider transitions between the ground state, 
I^m) ~ \-^) ® |1m)' where |1^) is a single-photon state, 
and the excited state, |e) — \X) ® |0), where \X) is 
given by Eq. ([5]) and |0) denotes the vacuum state. We 
must now solve the interaction picture Schrodinger equa- 
tion, j^\'li{t)) = ~-j^H'\^{t)), assuming that the system 
may be in a superposition of the two eigenstates, i.e., 
= Cc{t)\e) + T,^Cg^{t)\gf,) and by projecting the 
result onto either (e| or we obtain the set of equa- 
tions 



dt 



xe,,. yd3rV':,k.(r)A^(r)pVv,k.(r) (13) 



d 
dt 



E ^e'^^*x(kc,kv)ce(t) 



kc,kv 



xe^- / d^rCk (i-)^;(i-)pV'c,k.(r). (14) 



Let us first turn to the spatial integrals, which include 
the momentum operator p = — iftV. We assume that 



A^{r) and the plane- wave part of '0k (r) are slowly vary- 
ing on the length scale of the lattice constant, so that 
these functions can be evaluated at each lattice site r„ 
and taken outside the integral: 



d^i-V';k,(i-)AM(r)P^'v,k.(r) 



-,-ikc 



ikv-r"! 



— / d^r<ko(r)P'^v,kv(r) 

VUG JVC 



=6,, 



p,,J d-^re-*-^A^(r)e"'-- 



(15) 



where UC denotes integration over one unit cell with vol- 
ume Vuc and we have used the orthogonality of the Bloch 
functions. The last equation is obtained by noting that 
since the Bloch functions are periodic, the integral over 
UC is the same for all r„ and may be evaluated sepa- 
rately. Also, we have assumed that the Bloch functions 
depend only weakly on k, so that they can be evaluated 
at k = 0. The sum can then be converted back to an 
integral and finally we have defined the Bloch matrix el- 
ement as Pcv = /uc '^^'''^c,o(i")P^v,o(r) The Bloch 
matrix element is a material parameter, whose magni- 



tude evaluates ti 



,32 



moEj,{x) 



where E„ is the 



Kane energy^, which depends on the indium mole frac- 
tion X in the InxGai_xAs alloy. 

The summations over k- vectors can now be carried out. 
At this point it is advantageous to change notation to 
the electron-hole picture by substituting kc — >■ k^ and 
kv — > — kh. By insertion of Eq. ([T5|) and interchang- 
ing the order of integration and summation in Eqs. (|13p 
and (|14p the resulting equations take the form of inverse 
Fourier transforms and we have 



dt nmo ^-^ (jjfj. 



e^^-Pcv J d^rx*(ro,r,r)A^(r) 



dr hmo u)^ 



(16) 



(17) 



xe^-Pvc / d^rx(ro,r,r)A:(r). 



By integrating Eq. ([T7)) with respect to time, inserting 
the result in Eq. (jl6p . and finally rewriting the result by 
multiplication with a Dirac delta function in frequency 
and integrating over frequency we obtain 



■nq 



X / dt'e-'^-(*-*')ce(<'), 



(18) 



where the projected non-local interaction function is de- 
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fined as 



(19) 

where ep is tlic unit vector parallel to Pcv Equation 
[T8l is a main result of this work. It is valid beyond the 
Markov and DA approximations and in arbitrary optical 
environments; it is therefore a generalization of existing 
theories of dipole emitters in nanophotonic structures. 

If the term pNL(ro,'^)/'^ in Eq. is spectrally slowly 
varying over the linewidth of the emitter, we may eval- 
uate it at the emission frequency ujq and take it out- 
side the integral. In this Wigner-Weisskopf approxima- 
tion wc obtain -^Ce(i) ^ 



dae-"'^ = 2-k5{P) and da5{a) 



PNL(ro,c^o) ^^(^)^ ^^^^ 



i have been 

used. In the following we will write wq as w for brevity. 
By assuming that the exciton is initially excited (co (0) = 
1) we obtain the radiative decay of the exciton state pop- 
ulation |ce(t)P ~ 
rate is defined as 



-r(ro,w)t^ where the radiative decay 



r(ro,a.) 



(20) 



Since r(ro,a;) depends on the exciton envelope function 
through the projected non-local interaction function, it 
is not possible to state in general whether the decay rate 
will increase or decrease when calculated beyond the DA; 
the decay rate must be calculated for a given exciton 
state in a given dielectric environment. The physical sig- 
nificance of this result is clearer when expressed in terms 
of dyadic Green's tensors. This relation is derived in Ap- 
pendix [O and from Eqs. and (|C9p we obtain the 
important result 

PNL{ro,uj) |pcv|^y"d^r ^ d3r'x(ro, r, r)x* (ro, r', r') 



X (e^.Im{G(r,r',c^)}.ep) 



(21) 



The Green's tensor is a propagator of the electromagnetic 
field. Thus, G(r, r',aj) may be interpreted as the field 
amplitude evaluated at the position r due to a dipole at 
r' with frequency w. This is illustrated in Fig. [TJ To actu- 
ally calculate pNL(ro, w), we must obtain the Green's ten- 
sor describing the electromagnetic environment for the 
particular geometry. 

It is important to stress that the formalism developed 
here docs not change the selection rules for optical transi- 
tions, i.e., they are governed by the usual dipole selection 
rule according to which the change in angular momentum 
in the transition must be Am = ±1. This is fundamen- 
tally different from atomic quadrupolc transitions, where 
Am = ±2. The reason is that the quantum states of exci- 
tons consist of both an envelope and a Bloch part. Since 



the electromagnetic field at optical frequencies is slowly 
varying over a unit cell the approximation in Eq. (jlSp is 
very good. This unit cell DA could in principle break 
down for higher frequencies of the electromagnetic field 
and thus lead to multipole effects at the Bloch-function 
level but this is not relevant for the systems studied here. 
We note that a classical calculation of the dissipation rate 
of an extended dipole emitter leads to the same form of 
the non-local response as considered above; this is dis- 
cussed in further detail in appendix [B] 



IV. THE DIPOLE APPROXIMATION 

Before exploring the effects beyond the DA it is in- 
structive to consider the DA limit of the expressions de- 
rived above. In this case we can evaluate the Green's 
tensor at the center coordinate of the exciton, Tq. We 
obtain 

2 



PNL(ro,'-^) 



TTC 



r (e^ •Im{G(ro,ro,w)}- 
y d3rx(ro,r,r) 



(22) 



Obviously, the integrals over the envelope functions de- 
pend only on the excitonic degrees of freedom and it is 
therefore natural to redefine Eq. ([^0]) as 



DA 



(ro,^^) 



7rq 



hm?, 



0^0 



J d3rx(ro,r,r) 



/o(ro,w) 



(23) 

where we have introduced the familiar notion of the pro- 
jected LDOS^i-^, 



p(ro,w) 



2uj 



[el ■ Im{G(ro,ro,a;)} • 



(24) 



The LDOS is obtained by solving Maxwell's equations 
and it enters the quantum optical theory of light-matter 
interaction as the local density of vacuum modes that 
spontaneous emission can occur to. 

The interaction strength between an emitter and light 
can be characterized by the OS denoted /(w). We define 
this dimensionless quantity as the ratio of the radiative 
decay rate in a homogeneous medium rrad,hom(w) to the 
radiative decay rate rci(a;) of a classical harmonic oscil- 
lator of elementary charge^!, i.e.. 



where 



rci(w) 



r ad. horn 



(c) 
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67rmoeoC'^ 
We can rewrite the decay rate Eq. 



as 



rDA(ro, w) 



2TOQeo 



/(cj)p(ro,a;). 



(25) 



(26) 



(27) 
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where we have used Eqs. ([25|) and (|26p to obtain the OS 

2 



/(-) = 



y d3rx(ro,r,r), 



(28) 



which is independent of ro . The usefulness of the notion 
of the OS in the DA is apparent from Eq. ([?7)) . i.e., the 
decay rate is given by the product of the OS and the 
LDOS and thus the OS quantifies the strength with which 
the emitter interacts with hght. 

The OS can be calculated readily for the exciton mod- 
els discussed in appendix \^ In the strong confinement 
regime the result is 



/(-) = f 



J d3rFe(ro,r)Fh(ro,r) 



(29) 



where Fo(ro,r) and i^h(ro,i") describe the independent 
electron and hole envelope functions, respectively. Thus, 
we obtain the well-known strong-confinement resul^i^ii^ 
in which the OS is proportional to the overlap of the 
electron and hole wave functions. As opposed to the re- 
sult obtained in the single-particle picture^i^, there is 
no complex conjugation of either Fe(ro,r) or Jh(rojr), 
which is a result of the two-particle formalism used 
]-^gj.Q2j6|39 rpj^g wave function overlap integral in Eq. (23) 
cannot exceed unityi^ii^^ and hence the maximum OS in 
the strong confinement regime is given by /max = 

which shows that ^ can be interpreted as the OS of the 
bulk crystal without confinement and exciton effects. 

For the spherical exciton in the weak- confinement 
regime, 



(30) 



where L is the exciton radius and ao is the exciton Bohr 
radius. This is the giant-OS effect, i.e., the OS is pro- 
portional to the volume of the exciton. It is also strongly 
dependent on the exciton Bohr radius and therefore it 
shows a strong dependence on the effective masses of 
the carriers. For the disc-shaped exciton in the weak- 
confinement regime we have 



/(-) = sf^ 



(31) 



In this two-dimensional model the OS is proportional to 
the exciton area, which in the absence of inhomogeneities 
inside the QD^ is given by the area of the QD. 

In Fig. [1] we compare the calculated spontaneous emis- 
sion rate within and beyond the DA. Fig.[TjA) shows the 
DA result in which the decay rate is proportional to the 
LDOS, cf. Eq. 1^^. The classical interpretation of the 
LDOS is that it describes self- interference, i.e., it is the 
field strength at the center position of the emitter, ro, 
due to the emitted light. This is given by the propagator 
of the field, i.e., the Green's tensor G(ro,ro,a;), which 
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FIG. 2. Comparison between the radiative decay rates for 
spherical Ina;Gai_iAs QDs calculated within (dashed red 
lines) and beyond (solid blue lines) the DA for various emis- 
sion wavelengths A and indium mole fractions x as indicated 
in the figure. For large radii the decay rate is quenched. 



is indicated as an arrow in Fig. HJA). In the QED in- 
terpretation of spontaneous emission it is stimulated by 
vacuum fluctuations whose density is given by the LDOS. 
Spontaneous emission beyond the DA is governed by the 
double integral appearing in Eq. (pij) . i.e., it is given by 
the interference between all points in space weighted by 
the exciton envelope function as indicated in Fig. [TjB). 



At this point it behooves us to clarify the criterion for 
the validity of the DA. For a homogeneous medium the 
field distribution functions take the form of plane waves, 

ik-r 

i.e., A^(r) = ^^^^ , where e-c is the dielectric constant of 
the material and V is the quantization volume. From 
Eq. (IT9|) . it is clear that the DA holds when the field dis- 
tribution functions are slowly varying on the scale of the 
variations in the envelope functions. This is equivalent 
to the criterion |k|£ ^ 1, where C is the characteristic 
length scale of the emitter. In an inhomogeneous medium 
the field can be expanded in terms of plane waves. This 
means that there is not a unique k for which we can eval- 
uate this criterion. Thus, at a given frequency for which 
|k|£ <C 1 holds in a homogeneous medium, it will not 
hold in general for all k-componcnts of the plane-wave 
expansion in an inhomogeneous medium. This indicates 
that the use of the DA even for small QDs embedded 
in photonic materials, such as plasmonic nanostructures^ 
and photonic crystals, needs further justification. In gen- 
eral, one must simply compare the decay rate calculated 
in and beyond the DA to assess if it is valid. 
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SPONTANEOUS EMISSION DYNAMICS OF 
LARGE QUANTUMD DOTS IN SPECIFIC 
INHOMOGENEOUS MEDIA 



The non-local interaction function depends on a com- 
plex interplay between the specific geometry of both the 
electromagnetic environment and the exciton wave func- 
tion and physical insight into spontaneous emission be- 
yond the DA can be gained by considering the special 
cases discussed in this section. We calculate the radia- 

I 



five decay rate for spherical and disc-shaped excitons in 
homogeneous media in and beyond the DA as well as 
for disc-shaped excitons near a semiconductor-air inter- 
face. Beyond the DA the notion of the OS is less useful 
because light and matter degrees of freedom cannot be 
separated, i.e., Eq. ([27| is not valid beyond the DA. We 
could still use Eq. to obtain a dimensionless quantity 
characterizing the radiative decay rate in a homogeneous 
medium but the radiative decay rate in inhomogeneous 
media is neither proportional to the OS nor to the LDOS 
so we shall refrain from doing so. 

In the strong confinement model we obtain 



PNL(ro,w) = |Pcv|'^ J d3rFe(ro,r)Fh(ro,r) J dVi^,* (rg, r')Fh* (ro, r') (e^ • Im {G(r, r', 



^)} 



(32) 




X=75Q nm, x=0 



X=9Q0 nm, x=0.5 



H ^ h 



>.=1050 nm, x=1 



50 100 150 200 250 300 

Disc-shaped quantum dot radius L (nm) 

FIG. 3. Comparison between the radiative decay rates for 
disc-shaped In^^Gai-ajAs QDs calculated within (dashed red 
lines) and beyond (solid blue lines) the DA for various emis- 
sion wavelengths A and indium mole fractions x as indicated 
in the figure. In this case the decay rate saturates for large 
radii. 



In systems with pronounced anisotropy between electron 
and hole wave functions and large optical field gradients, 
this can give rise to a significant orientational dependence 
of the radiative decay rate even for small QDs, as was 
recently observed experimentally^. Here, however, we 
shall not explore this further because our focus is on large 
QDs. 

Let us now consider spherical QDs with parabolic con- 
finement potentials in an optically homogeneous medium. 
As shown in appendix [D] the non-local interaction func- 
tion can be evaluated analytically in this case and the 
resulting decay rate is 



r VTrEci^ — 



r 



Thus, the decay rate is proportional to the bulk crys- 
tal OS, the giant-OS term, (L/ao)^, and finally 
2c ) = g-(7rnL/A) ^ -yifhere A is the vacuum wave- 
length of the emitted light, which is an additional term 
originating from the breakdown of the DA. The compe- 
tition between these terms leads to a maximum in the 
radiative decay rate at 



where 



Vex 

27rn ' 



Ep ( 



(34) 



(35) 



(33) 



Equation ([33]) . is plotted in Fig. [2] along with the DA 
result, Eq. ([50)1 . In a microscopically realistic model the 
transition energy would depend on the geometry, size, 
chemical composition, and strain of the QDs. Here we are 
not concerned with such microscopic details and we sim- 
ply take the chemical composition and transition energy 
as being mutually independent and constant parameters. 
Thus, in Fig. [2] we vary both in realistic combinations as 
indicated in the figure. We assume that the refractive 
index of the QD can be approximated by that of the sur- 
rounding medium. We describe the surrounding medium 
as GaAs and include the frequency dependence of the re- 
fractive index as described in Ref. |40. For simplicity we 
consider only heavy-hole transitions and neglect the ef- 
fect of strain (the axial approximation) in which case the 
effective mass is isotropic, see Ref. [2O for further details. 
Figure [5] illustrates the results expressed by Eqs. ([M)) 
and (|35p . i.e., the radiative decay rate attains a maxi- 
mum when L ~ imax that depends strongly on ag and 
thereby on the indium mole fraction. This shows firstly 
that pure GaAs is a more promising material for achiev- 
ing a large light-matter coupling strength as compared to 
indium-rich alloys, due to the smaller exciton Bohr radius 
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of GaAs excitons, and secondly that for given material 
parameters, a fundamental limit to the light-matter in- 
teraction strength is imposed by the breakdown of the 
DA. A similar size-dependence has been predicted for 
ZnO QDs using semiclassical approaches —i^. Secondly 
it shows that even for small QD radii the DA leads to a 
systematic ovcrcstimation of the light-matter interaction 
strength. 

The fact that the radiative decay rate calculated be- 
yond the DA vanishes for a vanishing QD radius, cf. 
Fig. [5] is correct but here it arises for the wrong rea- 
sons. It is an artefact of the exciton model used here be- 
cause the weak confinement approximation breaks down 
for small radii, i.e., L 3> ctg is not fulfilled. For the param- 
eters in both Fig.[5]and Fig. [31 Cq attains a value of 12, 16, 
and 29 nm, respectively, in the subfigures (A), (B), and 
(C). In a more realistic confinement model but in the DA, 
the giant-OS effect reemerges for very small QD radii be- 
cause the envelope wave functions are strongly expelled 
from the QD, i.e., in this regime the excitons expand^ 
when the QD becomes smaller. From Fig. [2] we can thus 
predict that the QD size dependence of the radiative de- 
cay rate in a more realistic confinement potential would 
exhibit two maxima: one due to the giant-OS effect and 
its quenching for large QD radii as obtained in Fig. [2] and 
another giant- OS effect and its quenching at very small 
radii. For either vanishing or infinite QD radii the decay 
rate vanishes due to the breakdown of the DA. 

For the analysis in this paper we have assumed that 
the subband energy level spacing, A^, largely exceeds 
the thermal energy, ksT. For a relative effective heavy- 
hole mass of mhh = 0.59 and L = 300 nm the hole sub- 
band spacing is, cf. Eq. (jA6| . AEh — 62 mK so for the 
range of radii in Figs. [5] and [3] the experimentally required 
temperatures are accessible with standard dilution refrig- 
erators. We have also assumed that the level spacing ex- 
ceeds the homogeneous linewidth of the emitter, ftFiad- 
This criterion is not fulfilled for all values of L in Figs. [H 
[21 and m but could be valid in other materials in which, 
e.g., the transition energy would be higher. Beyond these 
approximations, several subbands would be populated^ 
and eventually the system would approach the bulk limit, 
which is beyond the scope of the present work. 

We have numerically calculated the radiative decay 
rate for disc-shaped QDs and the result is shown in Fig. [31 
along with the DA result Eq. (|3T|) . Here we consider a 
3nm thick QD, with varying lateral size. We use the 
same approximations as for the sphere considered above, 
except that here we include the anisotropy of the effec- 
tive mass relevant for a strained InGaAs layer embedded 
in GaAs as described in Ref. [2^ These results indicate a 
similar scaling of Lmax and /max as predicted by the an- 
alytical results obtained for spherical excitons. Since we 
keep the thickness constant, the OS does not vanish for 
large QD sizes as opposed to the sphere considered above. 
This is in agreement with results considering a non-local 
susceptibility of large quantum discsi^. Our calculation 
includes the numerical integrations also along the axial 




100 200 300 

Distance to interface z (nm) 

FIG. 4. Radiative decay rate of excitons confined in disc- 
shaped QDs embedded in GaAs as a function of distance to 
a GaAs-air interface for various QD radii L as indicated in 
the plot. For radii up to L = 32 nm the distance dependence 
resembles the LDOS because the DA is approximately valid 
but for larger radii the behavior changes completely. This is 
a striking effect of light-matter interaction beyond the DA. 



direction of the QD but the results are not changed sig- 
nificantly by assuming the DA in the z-direction. Also in 
this case the maximally achievable radiative decay rate 
for pure InAs QDs (Fig.[31[C)) is much inferior to that of 
GaAs QDs. 

Let us now finally consider disc-shaped QDs near 
semiconductor-air interfaces. The interface leads to re- 
flections, which alter the light-matter interaction. We use 
the same parameters as in Fig. ^A) and use the Green's 
tensor describing the proximity of the interface^Si^dl 
and the result is shown in Fig. [31 For radii up to 32 nm 
the oscillations coincide with the characteristic oscilla- 
tion of the LDOS^ apart from a small overall reduction 
in the decay rate, which is consistent with Fig. [31^ A), 
i.e., for small radii the DA overestimates the actual light- 
matter interaction strength. For L = 64 nm the oscilla- 
tion also appears similar to the LDOS, but the DA result 
(not shown) is about 50% higher than the result of the 
full theory carried out beyond the DA. For even larger 
radii (L = 128 nm and 256 nm) the decay rate oscillations 
change dramatically. In this regime, which is far beyond 
the validity of the DA, the spatial dependence of the ra- 
diative decay rate exhibits pronounced deviations from 
the LDOS and develops into a standing wave pattern. 

A comparison between Fig.[31and the highly non-trivial 
oscillations in Fig. [H leads to interesting implications for 
increasing the radiative decay rate. The giant-OS ef- 
fect provides an effective mechanism for increasing the 
radiative decay rate but this effect is quenched by the 
breakdown of the DA. By employing optical nanostruc- 
tures, the radiative decay rate can be enhanced beyond 
that limit. In fact, for the parameters of Fig. [31 A) the 
saturation occurs for a radiative decay rate slightly above 
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200 ns"'^ but near a semiconductor-air interface for the 
same parameters (Fig. the radiative decay rate can 
exceed 700 ns~ ^ . This is a direct example of the inter- 
twining of hght and matter degrees of freedom imposed 
by the breakdown of the DA. The semiconductor-air in- 
terface considered here leads to an increase in the radia- 
tive decay rate of more than a factor of three as compared 
to a homogeneous medium and exploring these effects in 
other nanophotonic structures such as photonic crystals 
or optical microcavities, where the effects could be much 
larger, would be a very interesting future direction of re- 
search. 



VI. CONCLUSION 



Appendix A: The exciton confinement regimes 

In this section three confinement regimes of excitons in 
nanostructures are discussed. The two-particle effective 
mass equation, Eq. ([S]) , cannot in general be solved ana- 
lytically for realistic QD geometries and heterostructure 
confinement potentials, but it can be solved within cer- 
tain approximations and more importantly in different 
limits of the ratio between the Coulomb energy and the 
conduction (valence) band subband spacing, Ai5o (Ai5h), 
in the absence of the Coulomb interaction. In the case 
where the Coulomb interaction is negligible, the exciton 
is said to be in the strong confinement regime and when 
it dominates over the confinement potentials, the exciton 
is said to be in the weak confinement regime^ii^iii^. 



We have derived the fundamental equations governing 
excitonic spontaneous emission beyond the DA. The DA 
cannot be assumed valid a priori in nanophotonic struc- 
tures even for small QDs and thus we calculated the re- 
sult for a solid state emitter beyond the DA. We derived 
the relation to the Green's tensor description of the elec- 
tromagnetic field. In this theory the radiative decay rate 
of excitons is described by a non-local interaction func- 
tion, which reduces to the LDOS in the dipole limit. The 
theory contains also the giant-OS effect in the weak con- 
finement regime as well as the strong confinement regime 
as limiting cases. We have investigated and clarified the 
conditions under which the DA is valid. We notice that 
these conditions depend on both the properties of the 
QD, the emission wavelength, as well as the structuring 
of the environment and the position of the QD. Thus, the 
DA is in general only valid in certain points of space for 
a given emission energy of the QD. 

Finally we note that the two cases discussed here, 
namely QDs embedded in either a homogenous medium 
or near a semiconductor-air interface, benefit from the 
availability of exact Green's tensors, simple experimen- 
tal realization, and therefore the possibility of direct com- 
parison between experiment and theory, but they are also 
the systems where the expected magnitude of the effects 
arising from the breakdown of the DA are smallest. A 
very interesting future direction would be to calculate the 
radiative properties of spatially extended excitons, e.g., 
in a photonic crystal cavity where very large differences 
between the DA and the theory developed here could 
arise. 



1. The unconfined regime 

In the absence of confinement, i.e., when ^^(ro) and 
Vh(rh) can be completely neglected, Eq. ([8|) reduces to 
the problem of a hydrogen aton>^ with effective masses. 
This describes a free exciton in a bulk semiconductor 
and in this model wc can calculate the characteristic en- 
ergy and length scales for an exciton. The characteristic 
length scale of the interparticle distance is given by the 
exciton Bohr radius. 



ao 



(Al) 



where eo denotes the vacuum permittivity, is the rela- 
tive static permittivity of the material, and the reduced 
mass is defined as 



momhh 

TO — . 

Too + TOhh 

The Coulomb potential is given by 
Vcoui(rc,rh) 



47reoer ro - rh 



(A2) 



(A3) 



so the energy scale, i.e., ionization energy associated with 
the free exciton is the effective Rydberg energy 



7^ = 



9 



47reoerao rnQma^ 



2. The strong confinement regime 



(A4) 
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In the opposite limit when Coulomb interaction may be 
neglected (strong confinement), the electron and hole are 
decoupled in Eq. ([8]), which reduces to two independent 
particle-in-a-box problems. These are readily solved and 
the solution is 

X(ro,ro,rh) = Fc(ro, ro)-Fh(ro, fh), (A5) 

where Fo(ro,ro) and Fh(ro,rh) are the electron and the 
hole envelope functions, respectively. Depending on the 
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purpose of the theoretical description, a simple model 
assuming infinite barriers, isotropic masses, and a simple 
geometry such as a cylinder or cube, may sufiice. In 
such cases an exact solution is readily available in the 
literature. If we consider a cubic QD with side length 
2L, the subband spacings are^^ 



2^2 



(A6) 
(A7) 



In the strong confinement regime, TZ <^ ae^+aEi, ^ ^j^j^j^ 

implies that L <^ \J ^f^ao ~ IAgq. This simple model 
with infinite barriers overestimates the barrier heights 
and we will in general use a heuristic definition of the 
criteria and assume strong confinement for L <^ gq and 
weak confinement for L ^ oq, where 2L is the spatial 
extent of the QD. 



3. The weak confinement regime: the spherical 
quantum dot 

For a spherical QD with parabolic radial confinement 
we have 



V;/h(ro,i-c/h) = X 



(A8) 



where the confinement potential is given by 51 and we 
define the radius of the QD as L = 2^ -j^^. We introduce 
the relative and ccnter-of-mass parameters^ 



R 

r 
P 

P 

M : 

m = 



mhrh 



rric + TTih 

Po + Ph 

"ihPc - n^oPh 



ruc + rail 
Tie + nih ' 



(A9) 

(AlO) 
(All) 

(A12) 

(A13) 
(A14) 



By this transformation the effective-mass Hamiltonian 
Eq. ^ separates into two decoupled Hamiltonians 



7JEM(ro,r,R) 
i?R(ro,R) 

Hr(r) 



FR(ro,R)+i/r(r) 
p2 ^ 



2moM 

v^2 



(A15) 

Mif\R~Yo\'^ (A16) 



(A17) 



2mom 47reoer|r| 

where we have neglected the term imfl^lrp in Eq. (|A17[) 
since we consider the weak confinement regime^. Thus, 
we have reduced the problem to solving the effec- 
tive mass equation for two well-known Hamiltonians, 
namely the three-dimensional isotropic harmonic oscil- 
lator, Eq. (|A16[) . and the hydrogen problem, Eq. (jA17p . 



We can then write the solution to the effective mass equa- 
tion as 

x(ro, To, rh) = x'(ro, r, R) = xcM(ro, R)Xrci(r), (A18) 

where XCM(ro7 R) is the center-of-mass wavefunction and 
Xrci(i") is the wavefunction describing the relative motion. 
For the present purposes, we are only concerned with the 
ground state envelope wave functions, which are given 

3/4 / , X 3/2 



XcM(ro,R) 



Xrcl(r) = 3 



TTttr 



1/2 



-kl/ao 



(A19) 



(A20) 



where P = \J and oq is the exciton Bohr radius. With 

these definitions we have L — \/2/? and by comparison to 
the definition of the normal distribution function we find 
that L equals two standard deviations, which we define 
as the radius of the QD. 



4. The wealt confinement regime: the disc-shaped 
quantum dot 

For a disc-shaped QD with harmonic in-plane confine- 
ment and infinite barriers in the z-direction we have in 
cylindrical coordinates, {p,z,4>), thatii^ 

K/h(ro,ro/h) = "^^Zc/h(^o,2c/h) + ^™c||/h||^^^|Po/h - Pop 



and 



VzSzG,z) = Vzjzo,^) 



(A21) 



Ofor|z-.o|<| (A22) 
oo for \z - zol > ^ ' 



(A23) 

where is the height of the QD. We assume that 
^ Oo so that the Coulomb interaction in the z- 
direction may be neglected. The infinite potential in 
the ^-direction is a somewhat crude approximation and 
in a more realistic model the wave functions would ex- 
tend into the barriers. However, if we model the system 
with Lz being slightly larger than the physical height 
it is a reasonable approximation although it does neglect 
the difference barrier penetration depths of electrons and 
holes due to the difference in their effective masses. This 
can be considered as a model of a quantum well with 
thickness fluctuation potentials^iilii^. 

For the in-plane coordinates, r = (p, </>) and p = 
(PptPip), respectively, we can make the same transforma- 
tions as in Eqs. (|A9|) to (|A14|> i^ with which the effective- 
mass Hamiltonian Eq. ([5]) separates into four decoupled 
Hamiltonians 

ifEM(ro,I", R, Zc^h) 

(A24) 
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where 



given bji 



16.24.45 



i^R(ro,R) 
i7r(r) 



2mom 47reoer|r| 



Pc 



2momc2 



-Mn^\Ii-ro\^ (A25) 

„2 

(A26) 
(A27) 



+ Vz,Szo,Zh), (A28) 



XcM(ro,R) 
Xroi(r) 

XzAzo,Zc) 

Xzi,{zo,Zh) 



2 1 

^g-2|r|/ao 

27rao 
2 / 7r(zc - 20) 

2 / 7r(zh - zo) 
■ cos 



(A30) 
(A31) 

(A32) 

(A33) 



Appendix B: Interaction between spatially extended 
classical emitters and classical electromagnetic fields 



TOf72|r|2 in Eq. (|X26)) 



where we have neglected the term ^ 
since we are considering the weak confinement regime^ 
Thus, we have reduced the problem to solving the ef- 
fective mass equation for three well-known Hamiltoni- 
ans, namely those of the two-dimensional isotropic har- 
monic oscillator, Eq. (jA25p . the two-dimensional hydro- 
gen atom, Eq. (|A26[) . and the particle in an infinite- 
potential box problem, Eqs. (|A27|) and (|A28|) . The solu- 
tion is 



x(ro,r,R, Zc^h) 

= XCM(ro, Tl)xi-c\{r)XzJzQ, Zc)Xzh(^o, Zh), 



(A29) 



where xcM(rojR-) is the center-of-mass wave function, 
Xvei(j) is the wave function describing the relative mo- 
tion, and Xzaizo,Zc) and Xzi,izo,z\i) describe the elec- 
tron and hole wave function in the z-direction, respec- 
tively. The ground state envelope wave functions are 

I 



It is instructive to consider spontaneous emission be- 
yond the DA in a classical model; this leads to an expres- 
sion, which is very similar to the quantum result derived 
in sec. IIIIl A classical emitter can be described by the 
current density 



J(ro,r) = -iw/xp(ro,r). 



(Bl) 



where fj, is the dipole moment and yo(ro,r) is the density 
of the emitter centered at Tq. This definition implies that 
the emitter is considered as a continuous distribution of 
infinitesimal dipoles. The power dissipation rate is 
given by Poynting's theoreni^i 



i j^d3rRe{J*(ro,r).E(ro,r)}, (B2) 



dW 

~dt' 2 jy 

where V denotes the volume occupied by the emitter and 
E(r) is the electric field, which is given in terms of the 
dyadic Green's tensor G(r, r',aj) as 

E(ro,r) =iw/iAio / d3r'G(r, r', w) • J(ro, r'), (B3) 
Jv 

where /Xq is the vacuum permeability. By combining 
these relations we readily obtain the normalized decay 
rate in an arbitrary dielectric environment as the ratio of 
power dissipation in the arbitrary structure to that in a 
homogeneous medium. The result is 



r(ro,Lj) ^ Jyd^rJyd^r'p{ro,r)p{ro,r')nl ■ Im{G(r, r',a;)} ■ 
ro(w) Jy d^r Jy d^r' p{ro,r)p{ro,r')nJ^ • Im {Go(r, r', w)} • n^' 



(B4) 



where denotes a unit vector in the direction of the 
polarization of the emitter and Go(r,r',w) denotes the 
Green's tensor in a homogeneous medium. The homoge- 
neous medium decay rate, To{uj) does not depend on rg 
because of the translational invariance of Grcens's tensor 



I 

in a homogeneous medium. In the DA we assume 

p(ro,r) = <5(r-ro), (B5) 



so Eq. (jB4p reduces to the simpler and well-known re- 
sult^i 



rDA(ro,a;) • Im {G(ro, ro, w)} • 



TooAi^) n?' • Im{Go(ro,ro,a;)} • n^' 



(B6) 
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Appendix C: Calculation of decay rates using 
Green's tensors 



where Xp(r,a;) are expansion coefficients. By combina- 
tion of Eqs. (fUT]) to (fU4|) we obtain 



Here we derive the relation between the Green's tensor 
and the vector potential. The Green's tensor G(r, r',a;) 
for the electric field is defined as^ 



VxVxG(r,r',a;) ^er(r)G(r, r', = I5(r-r'), (CI) 



where I is the identity matrix. Since the field distribu- 
tion functions describe the spatial part of the solutions 
to Maxwell's equations, they satisfy the wave equation 



V X V X (e^A^(r)) - -^e,{r) (e^A^(r)) = (C2) 



and the orthogonality relation 



J d3re,(r) (e^A^(r)) • {e^-A;,{r)) = S^^^,. 



(C3) 



Therefore we can expand the Green's tensor in terms of 
these functions 



^ ^ (e„A* (r)) ® (e^,AJr')) 



where (8) denotes the dyadic product. Using the identity 

(C6) 



lim Im , 



2uj, 



-{5{UJ - UJ^) - SiiO + LJ^)), 



multiplying with (e^^* (r)) (g) (e^^p(r')), and summing 
over all /i we obtain the useful relation 



(C7) 



Here we have discarded the unphysical delta function 
S{uj + ujf^). Now, from Eqs. (jCSP and ()C7P we obtain 



Im{G(r,r',w)} 



(C8) 



From this result the equivalence of Eq. ([T^ and Eq. (PT|) 
can be found directly by performing a series of op- 
erations on both sides of the equation. By project- 
ing onto Gp from left and right, multiplication with 
'PcvPx(ro,r,r)x*(ro,r',r')i and finally integration 



G(r,r ,0;) — N (6pj4^j(r )) , (C4) over both r and r' we obtain 



2lo 



IPcvl^y d^r J dVx(ro,r,r)x*(ro,r',r')(e^-Im{G(r,r',a;)}.ep) = 

r r (C9) 

|pev|'5^|e^-ep|^y d3rx(ro,r,r)A;(r) J dVx* (ro, r', r')A^(r')<5(c. - c.^). 



Thus, the right-hand sides of Eq. (fT9|) and Eq. (f2T|) are 
identical. 



Appendix D: Explicit evaluation of matrix elements 
for gaussian wave functions 

Here we show the explicit analytical evaluation of the 
non-local interaction function Eq. (|2ip for spherical ex- 
citons in the weak confinement regime. By insertion of 
Eqs. (|JT8)) . (Pl9)) . and (|X20)) in Eq. ^ we see that the 



integral to be solved is of the form 

J d'r J dV/(r)r(r')Im{e^G(r,r')-e„}, 

(Dl) 

where /(r) = /oe-l'-l'/'^' and = y/8Tr~^^^ao^ (3'^ . For 
a homogeneous medium, no generality is lost by choos- 
ing ro = and we have therefore suppressed Tq in the 
following. This allows for an explicit evaluation of the 
non-local decay function for homogeneous media. In the 
case of non-homogeneous media, we can always express 
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the Green's tensor, and hence the integral in Eq. (|D1[) . as 
the sum of a homogeneous part and a scattering part so 
this calculation is useful also for inhomogeneous media. 

Using the expression for the Green's tensor in homo- 
geneous medial^ we can rewrite Eq. ()D1[) as 



(D2) 



fc2 do? ) 47r 



/2 / d^r / A^v'e-^ e-^ 



(D3) 



in which ^ = |r— r'|, r = |r|, and k ~ |k| is the magnitude 
of the wave vector in the background material, and jo 
and /iQ denote the spherical Bessel and Hankel functions 
of the first kind, respectively. The spherical Bessel func- 
tion can be rewritten in terms of other spherical Bessel 
functions and spherical harmonics aa^ 



oo n 



n— m— — n 



(D4) 



in which the spherical harmonics are defined as 



47r(n + my. 



(D5) 



where P^(cos6') is the associated Legendre function. 
We now perform the angular integration over r' to find 

m ~ and n = 0. In addition, we note that -^-^ 

I — 7i • 
acts only on jo{kr). Following Ref. l50| we now write 



j„(fcr)y^(r) = r2™ in which case we may express the 
derivative in terms of raising and lowering operators 
and respectively, defined as 



and with the actions (for < \m\ < n) 



(D6) 



{n + m + 2){n + m+l) ^^^^ 



(2n+l)(2n-h3) 

(D7a) 



'(n-TO)(n-m - 1) +1 



4n2 - 1 



(n-m + 2)(n-m + 1) _i 
(27i + l)(2n + 3) 



{n + m){n + m. - 1) _^ 
4n2-l ""-1 



For the a — x term we find 



(D7b) 
(D7c) 

(D8) 



Only the D^D^ term results in non- vanishing terms after 
angular integration and we have 



(D9) 



= -l^" - -i^^i (DIO) 



The angular integral over vanishes, leaving only the 
rig term. In this way we obtain the final expression for 
I^^ as 



(Dll) 
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